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Abstract 

We consider cosmological particle production in 1+1 dimensional string theory. The process is 
described most efficiently in terms of anomalies, but we also discuss the explicit mode expan- 
sions. In matrix cosmology the usual vacuum ambiguity of quantum fields in time-dependent 
backgrounds is resolved by the underlying matrix model. This leads to a finite energy density 
for the "in" state which cancels the effect of anomalous particle production. 



1 Introduction 



The spontaneous production of quantum particles in curved backgrounds is a profound phe- 
nomenon which has inspired progress in fundamental physics for several decades. In particular, 
cosmological particle production is, according to inflationary cosmology, the origin of all observ- 
able structure in the universe. It is clearly important to understand particle production better 
in a context where the quantum nature of gravity is taken fully into account. The purpose of 
this paper is to attempt this in the case of the matrix model for two dimensional string theory. 

In the usual treatment of quantum fields in curved space-time cosmological particle pro- 
duction is the tangible consequence of an ambiguous vacuum: the vacuum state defined with 
respect to modes which are natural at early times typically contains particles when analyzed 
with respect to modes which are natural at late times (for a review see [1]). This is quite similar 
to Hawking radiation from black holes [2]. In string theory, holographic representations of grav- 
ity in terms of gauge theories have led to significant progress in our understanding of Hawking 
radiation [3, 4, 5, 6, 7, 8, 9, 10]. In the gauge theory description, there is a preferred time 
and so a preferred vacuum; and Hawking radiation results arises through the normal quantum 
mechanical decay of a highly degenerate initial state [4, 6, 11, 12]. It is natural to expect that 
string theory would provide a similar insight into the nature of cosmological particle produc- 
tion. This, however, has turned out to be rather difficult, particularly because of problems of 
formulating string theory in time-dependent backgrounds. 

Recent progress in two dimensional non-critical string theory has improved the situation 
somewhat [13, 14]. In this case there is a well understood holographic description — the matrix 
quantum mechanics of open strings — and also a closed string field theory, the two dimensional 
collective field theory of the eigenvalue density [15]. The holographic description has no space 
and has a unique time, while, in the closed string description, space arises from the space of 
eigenvalues [15, 16, 17]. Small fluctuations of the collective field then represent the perturbative 
states of the closed string theory (which in this case is a single massless scalar), while nontrivial 
time-dependent states, involving macroscopic numbers of decaying DO branes, correspond to 
cosmological evolution. Such time dependent states can be studied unambiguously using the 
matrix model description. 

The interpretation of time-dependent solutions in matrix theory as matrix cosmology was 
introduced recently by Karczmarek and Strominger [18] (the solutions themselves have been 
known for some time [19, 20, 21]). In the present paper we generalize their solutions by 
recognizing them as transforms of the ground state. This representation immediately allows 
the generation of several infinite families of cosmological solutions. Although our derivation is 
classical the solutions clearly exist at the quantum level. 

The fluctuations around the cosmological backgrounds are efficiently described by the col- 
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lective field theory of the matrix eigenvalues. The fluctuating field in fact reduces to a massless 
scalar field in two dimensions, a popular toy model for studying vacuum ambiguities. The 
novelty in the present context is that the vacuum state of the scalar field is inherited from 
the ultraviolet completion of the theory, i.e. the matrix model. This means the short-distance 
divergences of the energy-momentum tensor are cancelled, rather than subtracted as in usual 
quantum field theory. It also means a specific, finite, energy is associated with the static 
vacuum, namely the energy which, in the matrix model, is the standard one loop energy of 
the ground state. In matrix cosmology this ground state energy will contribute a novel time- 
dependent term in the EM-tensor which happens to precisely cancel the usual contribution 
from particle production. The VEV of the EM-tensor is thus identical in the initial and final 
states of the cosmology, a highly unusual situation. 

It is often confusing what the correct observables are in time-dependent string theory. This 
concern seems particularly acute in matrix cosmology where the universe tends to have dramatic 
initial and/or final conditions, such as the complete disappearance of space-time at early or late 
times. In the context of quantum field theory the problem with such non-adiabatic evolution is 
that the particle concept is not useful. In matrix cosmology there is a preferred vacuum, i.e. a 
notion of no particle state which is universal, and applicable even as space-time is disappearing. 
This enable us to discuss cosmological particle in a setting that would be difficult to analyze 
using ordinary quantum field theory. 

This paper is organized as follows. In section 2 we introduce the cosmological solutions 
to the matrix-model. In section 3 we first discuss the causal structure and the observables of 
matrix cosmology. We then construct the explicit modes of the fluctuating quantum fields and 
deduce the corresponding Bogulubov coefficients for particle production. Finally, in section 4, 
we compute the energy-momentum tensor of the model, using anomalies. 

Note added: as this paper was being prepared we received [22] which has some overlap with 
the present work, particularly section 3. However, the main points of the works are different. 

2 Matrix Cosmology 

In this section we introduce the matrix cosmologies of Karczmarek and Strominger [18], along 
with our generalizations. We discuss in turn the matrix model description, the collective field 
theory, and the fluctuations in collective field theory. 

2.1 The Fermion Phase Space Picture 

The holographic description of two dimensional closed string theory is in terms of singlet states 
of the quantum mechanics of a single N x N hermitian matrix M with some invariant potential 
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Tr V(M) which has a quadratic maximum. We will choose the value of V(M) at this maximum 
to be zero. The dynamics can be entirely recast in terms of the eigenvalues Xj(t) of the matrix 
M and interpreted as N fermions with positions Xi(t) in an external potential V(x). The 
double scaling limit consists of tuning the coupling constants involved in the potential and 
taking N — > oo so that the Fermi energy —fi F — > while the rescaled Fermi energy \x = (3N/j, F 
is held fixed. In this limit the coordinates and the momenta of the fermions may be rescaled 
such that the single particle Hamiltonian becomes the inverted harmonic oscillator 

The Fermi energy of this rescaled problem is Interpreting the double scaling limit as a 
continuum limit of string world-sheets one can identify the corresponding string coupling as 

9s = I- 

In the classical limit we can discuss the dynamics in terms of fermion trajectories in phase 
space (x,p). In the ground state, the Fermi surface is 

\p 2 - \* 2 = -V , (2) 

which implies that all states with p 2 < x 2 — 2/i are filled. In this paper we will concentrate on 
the region x < 0. As is known now, the excitations of the Fermi sea in the x > region simply 
represents a second massless scalar present in the OB theory [23]. 

The action corresponding to (1) has an infinite symmetry algebra, [24, 25, 26], whose 
generators are given by 

Wr a = \e^\x-p) r (x + p) s . (3) 

The Hamiltonian is one of these charges, h = —Wu. At the classical level the charges w rs satisfy 
the Poisson bracket algebra 

{w rs , w r > s >} PB = (rs' - sr')uv+r , -i,H-«'-i • ( 4 ) 

Our interest in as that it generates nontrivial solutions: starting with any classical solution, 
a new classical solution can be found by transforming with an element. In particular a 
static Fermi sea like (2) can be transformed into a time dependent Fermi sea by using charges 
with r s. We will concentrate on solutions generated in this way by charges uvo an d Wq s . 
For these the finite transformations of x and p become 

w 0s : x' = x + \se~ st (x + p) 3 ' 1 , p' = p - \se~ st (x + p) s_1 . (5) 
w r0 : x' = x + \re rt (x — pY" 1 , p' = V + \re rt (x — p) r_1 , (6) 
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where A is the parameter of transformation. Thus, starting from the static Fermi surface (2), 
we can obtain an infinite set of exact solutions characterized by Fermi surfaces 

w 0s : ^{x 2 -p 2 ) + Xse- st {x + p) s = \i . (7) 

w r0 : -{x 2 -p 2 ) + Xre +rt {x - p) r = fi . (8) 
2 

Of particular interest are the solutions obtained by the actions of woi and wio. In these 
cases the transformations (5-6) reduce to conventional coordinate transformations which 
are simply time dependent shifts of the coordinate. It is then clear that we can create a two- 
parameter family of solutions by combining the transformations Woi and w±o with arbitrary 
parameters 

x' = x + A_e~* + A + e* , (9) 
p = p- A_e~* + A+e* . (10) 

The resulting solutions 

l -{x 2 - p 2 ) + \_e~\x + p) + \ + e\x -p)=n, (11) 

are in fact the solutions discussed by Karczmarek and Strominger [18] (up to a redefinition of 
the string coupling constant). 

The more general cosmological solutions given by (8) also correspond to smooth Fermi 
surfaces. A nice way to map them is to introduce the coordinates x± = x ± p in phase space. 
Then the solutions generated by wo s take the form 

x.x + + 2\se- st x s + = 2/i . (12) 

We will consider the case of A > 0. After the further rescalings y_ = X-/2fia and y + = ax + 
with the time-dependent factor 

a = e -< , (13) 

the Fermi surface becomes 

y- = ^-(i - y s + ) ■ (14) 
y+ 

Interestingly this surface is symmetric around y = in phase space only for even s, since only 
then y__(—y + ) = —y_(y + ) . Concentrating as usual on the region in which y + < we find, for 
odd s a smooth Fermi surface interpolating between y^ ~ for small y + and y_ ~ y^ 1 

for large y + . Importantly, since the quantity vanishes only once in this domain, there are 
no "folds" in the Fermi surface. This means it can be characterized by its intersection with a 
x + = constant line, i.e. by the function a:_ = x-(x+). For even s, there is no zero of ^ and 



the Fermi surfaces actually cross over to positive values of y_. In this description the Euler 
equations for the Fermi surface are simply 

d t x± = ±x± ± x tq^T > ( 15 ) 

which may be easily verified for (7-8). Alternatively, we can parametrize the Fermi surface by 
the values of p at the two intersections with the x = constant line, P±(x,t). Then the Euler 
equations take the form 

d t p ± = x - p±o x p ± . (16) 

In the remainder of this paper we will primarily discuss the original solutions (11) of [18]. 
2.2 Collective Field Theory 

So far our discussion has been in the open string language of matrix quantum mechanics. From 
this point of view we have described some matrix configurations which are time dependent, but 
there has not been a notion of "space", and therefore no "cosmology" to discuss. The spatial 
coordinate x is an emergent quantity which can be seen only after passing to the closed string 
description by introducing the density of eigenvalues, or the collective field 

p(x,t) = Y,6(x-Xi(t)) ■ (17) 

i 

In the continuum limit it is convenient to trade the density p(x, t) for a scalar field <f)(x, t) 
through p(x,t) = d x (f)(x,t). The action of the field <j>(x,t) is [27, 28] 



S — J dtdx 
where AS is the singular term 



2d x (f) 6 v ; v 2 p; 



+ AS, (18) 



AS = i J dtdx d x (p 



d x d x i log \x — x'\ 



(19) 



which is a part of the Jacobian of the change of variables (17) from eigenvalues Xi(i) to the 
density p(x,t) [27]. The singular term AS contributes at higher order in the loop expansion 
parameter g s = p~ l , but it will nevertheless play a central role in what follows. 

Any distribution of eigenvalues corresponds to a definite state of the two dimensional col- 
lective field theory. However, for generic Fermi surfaces such a state cannot be described as a 
classical solution of the collective field theory because of the presence of folds or disconnected 
pieces [29, 30]. In fact, it can represent states of the theory where quantum dispersions of 
fields are of the same order as their classical expectation values [31, 32]. However, for the time 
dependent Fermi surfaces considered in this paper the profiles are quadratic — there are no 



folds — and such Fermi surfaces can indeed be represented as classical solutions of collective 
field theory. 

The classical equation of motion following from the Lagrangean (18) is 

2 



28,^ = a, 

o x d 



^ + m) - ,ji - 2 "» 



(20) 



since we must ignore AS in the classical limit. The ansatz 



d x h = -P (x, t) , d t (j> Q = —P (x, t)F(t) , (21) 

71 71 

where 

Po(x,t) = - F) 2 - 2/i , (22) 

solves the equation of motion (20) for all functions F(t). Imposing the consistency conditions 
d t d x (po = d x d t (fio on (21) we find 

and so 

F(t) = A_e-* - A+e* . (24) 

An alternative procedure that yields the solution (21) with (22) and (24) is to verify that the 
profiles of the Fermi surface 

P ± (x,t)=±P (x,t) + F(t) , (25) 

satisfy the Euler equations (16). 

The two parameter family of solutions to (23) given in (24) is identical to the matrix model 
solution (11) that was generated from the static solution F = by the action of w 01 and w w . 
The explicit solutions for other w 0s , w r0 solutions can be also obtained in principle, though they 
involve solutions of higher order algebraic equations. 

The form of (22) restricts the x < branch of interest to 



x < -y/2fi + F . (26) 

The solution generated by w i has F(t) = A_e~* so, as t = — oo, the Fermi surface is pushed 
to the region of large negative x. In perturbation theory, the collective field therefore does not 
support excitations for any finite x — there is no universe at t — — oo. At finite times there is 
some allowed region of x which is growing at t increases so that, as t — > oo, the Fermi surface 
becomes the static solution (2), which is the usual universe of the two dimensional string around 
the ground state. The entire solution thus represents the creation of a universe. In a similar 
way, solutions generated by w w represent the destruction of a universe. We will henceforth 
concentrate on universe destruction. 
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It will be important for our considerations to introduce a lower bound x niin = —A in (26) 
which restricts the universe to the finite volume 

-A < x < -J^i + F . (27) 

The origin of the infrared regulator is the potential V(x) before performing the double scaling 
limit. This leads to A ~ \fj3N ~ l/y/flo — > oo. We should therefore think of the volume of the 
universe as a quantity of order \f/3N in x space. 

The solutions of collective field theory which correspond to the more general solutions gen- 
erated by w 0s or w r0 are more complicated to obtain, since these correspond to nonlinear trans- 
formations of the collective field. Nevertheless, as mentioned above, they can be obtained by 
solving higher order equations in terms of phase space variables and translating these solutions 
in the collective field language. 

2.3 Backgrounds in Matrix model 

What is the meaning of these cosmologies in the original holographic theory, i. e. matrix quan- 
tum mechanics? This is an interesting question because in this case the holographic theory has 
no "space" at all and these cosmologies are encoded as specific time dependent configurations 
of the matrix M(t). The Hamiltonian for matrix quantum mechanics reads 

H M = Tr[U 2 M - M 2 ] , (28) 

and the charges are given by 

w rs = V~^Tr [(M - U M ) r (M + U M ) S ] . (29) 

If denotes the ground state of this system, a cosmological background generated by w s o is 
denoted by 

| A) = exp{i\w s0 ]\fi) , (30) 
The expectation values of a typical invariant quantity may then be formally expressed as 

(A|0|A> = (/i|OV) , (31) 

where 

qi = e -i\w s o 0e i\w s0 _ ( 32 ) 

This computation is purely formal because the left hand side of eq. (31) is not meaningful 
because the state |A) is not normalizable. States such as |A) form a subset of the discrete states 
in the c = 1 model [33]. These, it is believed, should not be really be regarded as states in the 
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spectrum of the original model. Instead, they are interpreted as deformed backgrounds. The 
formal manipulation suggests that the appropriate deformed Hamiltonian is 

H' = H M + Xse st Ti{M - U M ) S ■ (33) 

In other words, the holographic interpretation of the cosmological background is the modified 
matrix model whose Hamiltonian is H' . This interpretation is similar to the proposed descrip- 
tion of two dimensional black holes in terms a matrix model deformed by a TrM -2 potential 
[34]. 

2.4 Fluctuations 

The sole perturbative closed string excitation around an arbitrary classical background is given 
by the fluctuation of the collective field around the corresponding classical solution of collective 
field theory <p (x,t) 

t) = Mx, *) + ~^v(x, t) . (34) 



In our applications we identify the field rj with the spacetime tachyon 1 . Inserting the expression 
(34) into the action (18) we find the quadratic action of the fluctuations 



sm = ll 



dxdt 



(d tV f - LdM 2 - (fx) 2 ) (^) 2 - 2 inr d ^ 

V Ox<P0 J O x (po 



- , • (35) 

This form of the quadratic action is valid for fluctuations around any classical solution <j>o(x, t). 
It is useful to interpret the action in terms of a massless scalar field propagating in a spacetime 
with metric 

\2 



ds = ~ dt + f T^fk ■ (36) 
For the specific case of matrix cosmology (21) we have 

S W = I f dxdt — I— Ud t r,) 2 - (P (x, t) 2 - F{t) 2 ) (d xV ) 2 + 2F(t)d tV d x r ] \ , (37) 

and the background metric becomes 

The coordinate x takes values in the interval (26). At the end-points of this space the field 
f)(x,t) satisfies Dirichlet boundary conditions. This follows from the fact that the integral of 
the full collective field / dxd x (f) is the total number of fermions and therefore fixed 2 . 

1 In general leg-pole factors must be taken into account. These result in a nonlocal redefinition of the field 
which will not be essential in our discussion. 

2 Strictly this argument implies only that ?y(x m i n ) — r/(— v / 2/I)=0. That fermions are in fact prevented from 
leaking out in either end of the interval is clear prior to taking the double scaling limit, and this property is 
inherited by the scaled theory. 
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At the classical level the action (37) is invariant under Weyl rescalings and so it determines 
the metric only up to an overall conformal factor. In (38) this factor was chosen so that 
spacetime is flat. This can be seen explicitly by transforming the spatial coordinate as 

x = — \J2/1 cosh y + F , r — t, (39) 

so that the metric becomes 

ds 2 = -dt 2 + dy 2 , (40) 

and the boundary condition simply becomes a Dirichlet condition at y — 0. In these coordiantes 
the problem thus returns to the static case. 

In terms of the original coordinates (x, t) we have a problem similar to the moving mirror 
problem [1, 35], with the mirror trajectory given by the upper limit of the coordinate range 
(26). In the moving mirror problem particle production can be entirely rephrased in terms of 
the anomalous transformation of the energy momentum tensor under the conformal transfor- 
mation which makes the mirror stationary. The vacuum in static coordinates then appears as 
a collection of particles in the frame where the mirror moves. Particle production in matrix 
cosmology works similarly except that, as we shall see in section 4.2, the vacuum is prescribed 
differently. 

In usual quantum field theory there is, in the absence of special symmetries, no sense 
in which one vacuum is preferred over the other — the appropriate choice is determined by 
the nature of observers. This vacuum ambiguity enters concrete computations through the 
normal ordering prescription which is implemented at "equal time", a notion that depends on 
the observer or, more precisely, on the coordinate system. String theory is different because 
ultraviolet divergences are absent and therefore physical observables, such as the one loop free 
energy, are finite. The way this comes about in the collective field theory description is that the 
singular term (19) acts as a counterterm that cancels all divergences. The form of this singular 
term singles out a specific coordinate system, by regulating physical quantities at equal matrix 
time. This circumstance forces us to use the same time t — r in the original (x, t) coordinates 
and the static coordinates (y,r). 

3 Particle Production in Matrix Cosmology 

In this section we discuss the causal structure of the metric seen by the fluctuations and define 
notions of in and out modes accordingly. We derive the nontrivial Bogolubov transform relating 
the in and out modes, and we determine the corresponding spectrum of particles. We consider 
for definiteness the destruction of a universe, i.e. a draining Fermi sea (F = — Ae*). 
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3.1 Causal Structure 



A general problem in the study of time-dependent backgrounds is the definition of proper 
physical observables. For example, it is often difficult to define an S-matrix, because there are 
no suitable in and out regions. One might expect that matrix cosmology would suffer from this 
problem, since the entire spacetime disappears at late times. However, as we discuss now, this 
is fortunately not the case. 




A i B 

Figure 1: Causal structures in (A) static (y) and (B) cosmological (Q) coordinates. The dashed 
lines refer to the boundary of spacetime (the endpoint of the eigenvalue distribution) in the 
two different coordinates. 



The most natural coordinate for spacetime processes is not the matrix model coordinate x, 
but rather the exponentiated coordinate Q defined through 3 



x = -y2/icoshQ (41) 

We will think of Q as the physical spacetime coordinate and refer to it as the cosmological co- 
ordinate. In our applications we will in addition find it essential to employ the static coordinate 
y which transforms the matrix cosmology to Minkowski space (40) with a Dirichlet condition 
at y — 0. According to (39) we have the relations 

F A 

cosh y = cosh Q -\ — -== = cosh Q 7S= e * (42) 

y2/i v^/ 1 



3 In the literature one often encounters the coordinate q defined through x = — e q (see e.g. [39]). Our Q ~ q 
at asymptotic distances but Q has the advantage that the turning point of the Fermi sea is at Q = 0. 
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The time coordinate t is the same for all these coordinate systems. 

In the static coordinates (y) a natural in region X~ is defined by taking the light-cone 
coordinate y_ = y — t — > oo. Similarly a natural out region X + is defined by taking the light- 
cone coordinate y + = y + t — > oo. It is important to distinguish these regions from the future 
and past time-like infinities, i ± , which are reached by taking t — > ±00 at fixed y, and also from 
the space-like infinity i°, which is reached by taking y — > +00 at fixed time t. These are all 
indicated in figure 1(A). 

This discussion of asymptotic regions in the flat coordinates translates nicely into the cos- 
mological coordinates (Q). Indeed, define light-cone coordinates Q± — Q ±t and use (42) to 
obtain 

Q± = Q±t 

= ±t + costT 1 < coshy + 




± H±^ + cosh- cosh («±±±) + -A=e^ ,43) 



The in region X is defined by taking y_ — > 00 with y + fixed so (43) gives 

Q±=y± + 0(e' y -) on J" (44) 

The in region X~ is thus defined in physical coordinates by taking Q_ — > 00 with Q + fixed. 
Of course this is rather obvious, since it is clear from (42) that the static coordinates (y) agree 
with the cosmological coordinates (Q) at early times. The point of making the notion of "early 
times" precise is to avoid confusion about the concept "late times" in the following. 
The out region X + is defined by taking y + — > 00 with y_ fixed; and so (43) gives 

g ± = y± + log(l + Ae- J/ -)+0(e- 2/ +) onJ+ (45) 

where A = y^A. Thus Q + — > 00 with Q_ fixed will take one to the out region X + . The relation 
(45) is easily inverted on X + to find 

y ± = Q ± + \ og (l-\e-Q-)+0(e- Q +) (46) 



An important feature of this expression is the branch cut at Q- = log A. Thus we are reaching 
the out region by taking Q + — > 00 with Q_ fixed at some value larger than the minimal value 
(Q-)min = log A. The significance of this limiting value is that the boundary of the draining 
Fermi sea asymptotically approaches the light-like trajectory parametrized by Q_ = (Q-)min- 
This is indicated on Figure 1(B). 

In summary, the disappearance of spacetime by definition prevents the presence of a future 
timelike infinite i + . However, there is still a well-defined out region X + , defined by taking 
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Q + — > oo with Q_ fixed at some value larger than the asymptotic trajectory of the tachyon 
wall (<3-)min = l°g A. This means it makes sense to discuss observables in the form of S-matrix 
elements. 

3.2 The in and out Modes 

The equation of motion for the fluctuating collective field 77 can be solved exactly. This is most 
easily done in the static coordinates (y) where the equation of motion is satisfied by simple 
plane waves. After taking the Dirichlet boundary condition at y = into account, an obvious 
basis of solutions is given by 

v>w(y, t) = -±=e** sin uy (uj > 0) (47) 



and their complex conjugates. The normalization of (47) has been chosen such that these modes 
form an orthonormal basis with respect to the Klein-Gordon norm 



\Uu,> 



,u u )=i / dW (uo^u*, - u*#9 M « w ) =5{uj- uj') (48) 



on any Cauchy surface E. We can find the exact modes in the physical coordinates (Q) by 
solving (42) for y in terms of Q and t, and substituting the result into (47). The modes obtained 
in this way will in general have a complicated dependence on the time t, but this dependence 
become quite simple in the in and out regions. 

In the in region (X~) the static coordinates (y) coincide with the cosmological coordinates 
(Q) and so we can write (47) as 

u™(Q, t) = -Le- JQji sincuQ = ^ =e ^ Q n T~ (49) 

In the second equality we omitted the term depending on Q~ because, on X~ , this term does 
not contribute to the probability current. Equivalently, the inner product (48) reduces on Z~ 
to 

(ua/,u w ) = i J d Q+ (uu,d Q+ ul, - u^9q + u w ) on J" (50) 

According to (49) the modes (47) reduce to the standard positive frequency plane waves in the 
in region. It is for this reason that we have specified them from the outset by the superscript 
"in". 

In the out region (J + ) the static coordinate (y) is related to the cosmological coordinate 
(Q) through (46) and so the modes (47) take the form 

t£(Q,f) = -j=f**- (l " Ae-«-) i£J on J + (51) 



TILO 
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We have omitted the term depending on Q + because this term does not contribute to the 
current in the out region 4 . The canonical modes (49) in the in region thus evolve to the more 
complicated modes (51) in the out region. The change from dependence on Q + to dependence 
on Q- is due to the reflecting boundary conditions on the field which turn left-movers into right 
movers. 

Since the in modes are rather complicated in the out region it is natural to introduce 
a different basis which is simple there. The obvious choice is to consider a set of modes 
parametrized by uj > and which, in the out region, reduce to the canonical form 

u°J t (Q,t) = --^e^- onJ + (54) 

It is obvious that modes in fact exist that satisfy the equations of motion everywhere and reduce 
to this expression in the out region: the m-modes and their complex conjugates all satisfy the 
equations of motion and, in the out region, they take the form given in (51) which, when the 
complex conjugates are included, span all functions of <5_. The nontrivial content of selecting 
the modes (54) is the implied notion of positive frequency uj > which, after quantization, 
amounts to the introduction of a particle concept. Since spacetime is rapidly evolving at late 
times it is not obvious a priori that any such notion should even exist in the out region. However, 
in the present context, the underlying matrix model singles out a preferred time coordinate t. 
The proposed out modes (54) are uniquely determined by their canonical dependence on this 
time coordinate. 

3.3 The Bogolubov Transformation 

The in modes depend on both the out modes and their complex conjugates, as encoded in the 
Bogolubov transform, 



= I *> H^K"'(Q-) + «"'«-)) ] • (55) 

We want to compute the Bogolubov coefficients a(uo,uo) and (3(lj,uj). To do so we introduce 
the Fourier transform F of the in mode through, 

/oo 
dJ e iujQ -F(u,J) . (56) 
-oo 



4 This is clearest in the y-coordinates where the inner product on I + is written as in (50) with dQ+ — > dy- 
and 8q + — ► d y _ . If we wish to write the inner product on I + as an integral over dQ- we must use the more 
complicated tangent derivative 



d 



0Q+ 



dQ+ \e-Q- 



d Q+ (52) 

v+ 



dQ 

where 

(53) 



dQ- y+ 1 - \ e Q- 
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We must be a little careful when inverting this expansion, because has a semi-infinite range 
bounded below by the minimum value (<5-)min = In A. Accordingly we compute, 



f 

Jin 



A 27T 



In A 2lX J-oo 

00 dJ F(uj,uj' 



= 1 



(u>'—u>) In A 



2n uj' — uj 



= F(uj,lj) . 



(57) 



In the first step we assumed lm(u) < 0, which is implied already on the left hand side to ensure 
convergence. In the next step we closed the w'-contour in the lower part of the complex plane. 
Using (51) for the in-modes we find 



F(u,u) 



dQ. 



-iu>Q- ..in 



In A 27T 

% roc 



= / 

AlTUJ Jin 



\J Atiuj Jhi\ 2ix 



V Atxuj 2n 



Jo 



dz z -i( u -G>)-l (1 _ 



V 'Atiuj 2tt 



B(i(uj — cu), 1 + iuj) 



(58) 



We changed the integration variable to z — Ae _< ^~ . Comparing the definition of the Bogolubov 
transform (55) with the Fourier transform (56), and recalling the definition of out modes (54), 
we now find the Bogolubov coeffcients, 



a(u,u) = ^J-X^-VBMQ-w^l + iuj), 

2tT y UJ 

P(uj,uj) = ^\-X i{w+a)) B(-i(uj + uj),l + iuj). 
2ix V uj 



(59) 
(60) 



The non- vanishing of the (3{uj,uj) is interpreted as particle production in the matrix cos- 
mology. The creation and annihilation operators in the in and out vacua are defined by the 
expansions, 



77 ( 



duj 



r 

i(Q,t) = / 
J 

= duj 
Jo 

and so, comparing with (55), one finds, 



a™u™(Q,t)+a^u™*(Q,t) 
a: ut u°f(Q,t) + a^u°f*(Q,t) 



^out 



f 

Jo 



dw 



a(uj,w)a™ + (3*(uj,w)a, 



in\ 
w 



(61) 



(62) 
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The quantity that we want to compute is the expectation value of the out number operator in 
the in vacuum, 



N(w) = (OOm 

roc 

|2 



out] out\ 

dw \fl(w,w)\' 
00 dw sinh7ru; 
lo 4:7r(cu + w) sinh -k{oj + w) sinh ww 

This is our final result for the particle spectrum on X + . The integral is convergent at high 
energies, as it should be. There is a logarithmic divergence for small 00 which is interpreted as 
a linear divergence in Q_, i.e. it is the energy density which is finite. 

4 The Energy Momentum Tensor 

In this section we analyze the energy-momentum tensor of the fluctuations in the cosmological 
background. We begin by reviewing the static case and then turn to the time- dependent setting. 

4.1 The Static Case 

We begin the discussion by considering the static case F — 0. The Hamiltonian for the fluctu- 
ations then simplifies to 

H = J dxP ]±n 2 v + i(d^) 2 ] + AH, (64) 

where 

AH = J dxP 



— d x d x i log \ x — x'\ 
2ir 



(65) 



is the Hamiltonian form of (19). This last term is actually independent of the fluctuating field 
77, but it is included because it contributes at the same order in the loop expansion parameter 
g s = fi" 1 as the quadratic fluctuations. The measure dxPo that appear in both (64) and (65) 
translates in the underlying matrix model to dep(e), which is the obvious continuum version 
of the sum over eigenvalues. The expressions in the square brackets then correspond to the 
quantum fluctuations of the eigenvalues. 

For explicit computations it is useful to employ the static coordinates y where the metric 
for the fluctuation of the collective field is the simple Minkowskian (40), and the Dirichlet 
condition is imposed at y — 0. In this coordinate system the Green's function takes the simple 
form 

4vr & \At 2 -(y + yf 

where At = i—t and Ay = y — y. The denominator arises from the image charge that enforces 
the Dirichlet condition. The Green's function in other coordinate systems can be found by 
simply substituting expressions for r and y as functions of those other coordinates. 
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D(i, y- 1, y) = -— log ( A , 2 * , (66) 



The Green's function diverges at coincident points so, to do the calculations, we must reg- 
ulate the theory. We interpret two point functions by using a cutoff in x coordinates according 
to the prescription 

((d xV ) 2 ) = \imd x d x D(t,y;t,y), (67) 

where the point-splitting x — > x + |; x x — | is implied. The t and y (and their barred 
analogues) are functions of these slightly shifted x (and x). Evaluating the remaining expression 
we find 

_ , 2 \ 1 ,. f y'y' y'y' 1 

o x rj) ) = -— hm 



2 



27T [ (y - yf (y + y) 

' ~ ( ' A ~ 1 {//••'•} • O(r). (68) 



2vre 2 8tt \ y J 12tt 
where the limit e — * is implied and we introduced the Schwarzian derivative 



<*->-?-§($)'• 

as well as the notations y' = d x y = —1/Po etc. (and similarly for y and x). 

The canonical momentum is given by the Hamiltonian equation of motion as 

n„ = ^d t rj. (70) 

-M) 

Using this expression, it is straightforward to compute the two point correlator of the momenta 
using the Green's function (66) and the point-splitting procedure already used in (67). The 
result is 

This result agrees precisely with (68), except for the sign of the second term, which is the 
part that arises from the images. Indeed, the computations leading to the two results are 
almost identical because, after taking the Pq appearing in the denominator of (70) into account, 
temporal derivatives act on the short distance part of the Green's function in the same way as 
spatial derivatives. This agreement between the potential energy (68) and the kinetic energy 
(71) is referred to as the virial theorem. The origin of the virial theorem in the matrix theory 
is the simple oscillator form of the potential after double scaling. 

Collecting the results we find that the divergent pieces cancel and the final result for the 
expectation value is 

H gs = -J dxP -±-{y,x}. (72) 
of the Hamiltonian (64) in the ground state. The explicit form (39) (with F — 0) of y(x) gives 

{y ' x}= i? ^ + ^l)^2^' (73) 
16 



for large x 2 . The extensive part of the energy is therefore 

Hgs = ~i l ln (-^min)| = ^ ln ^o, (74) 
| = A ~ 1/ ^ S /J^Q. This result is interesting for several reasons: 

1. The singular term from the Jacobian (65) cancels the regularized singularities from the 
two-point functions and so the final result is finite, without the need for futher renor- 
malization of the collective field theory 5 . This explicit understanding of how the matrix 
theory induces the correct counterterms in collective effective field theory is the origin of 
a preferred vacuum in the theory. 

2. The image in the Green's function (66) contributes to the finite part of the two-point 
functions (68) and (71), but these contributions cancel in the total energy. It is not 
surprising that boundary conditions are unimportant for the extensive part of the energy 
in the thermodynamic limit; but it is nice to see how it works explicitly. 

3. The finite part agrees with the one-loop result found by solving the matrix model explicitly 
[36] 6 . This gives great confidence that we have interpreted the theory correctly. 

Let us also mention the further results: 

1. In collective field theory, the ground state of the matrix model is the vacuum defined in 
terms of the mode expansion of the fluctuation field in terms of modes e ±lbJt sm(uy) since 
this is the vacuum which leads to the two point function used above. 

2. The presence of a finite ground state energy in this theory is an important signature that 
we are dealing with a string theory. In fact, at finite temperature T this term, when 
added to the standard thermodynamic contribution, leads to a T-dual answer symmetric 
under ixT — > ^ characteristic of string theory [15]. 

3. We have performed the calculation using a point splitting regulator in x space. If we use 
instead a regulator in y space, the Schwarzian derivative term will come from AH. The 
divergences cancel as before leading to the same finite answer. 

5 This cancellation has been known from early days of collective field theory [15, 28, 37]. However, generally 
the image term has been ignored in these computations and our explicit verification of the virial theorem in the 
presence of regularization is also new, to the best of our knowledge. 

6 The energy (74) agrees with eq. 3.30 of [38] after dividing the result given there by two, because we only 
compute the energy of the x < part of the Fermi sea. 
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4.2 The EM-tensor in Matrix Cosmology 

Let us now generalize these considerations to matrix cosmology. As we have emphasized, matrix 
cosmology is related to the static case by a simple coordinate transformation which, near the 
asymptotic null infinity X + , takes the form 



y- 



= Q_ + log[l-Ae- Q -] + 0(e- Q +), (75) 



as Q + — > oo with Q_ fixed at some Q_ > In A. As in previous sections we use the notation 
y± = y ±t and Q± = Q ± t for the static and cosmological light-cone coordinates and we also 
introduced the abbreviation A = Ay^ to parametrize the limiting value (<5-)mm = In A. 

The static coordinates (y) play a several roles in the problem. First, the cosmological 
coordinates (Q) reduce to the static coordinates at early times; so the static coordinates define 
the in vacuum. Second, the Green's function (66) in static coordinates transform simply into 
the cosmological coordinates, since the transformation (75) is conformal (it does not mix the 
" + " and " — " light-cone coordinates). This will allow us to use conformal techniques to 
study the problem at asymptotic null infinity X + . Following the standard strategy, we begin 
the discussion by writing the results from the static case in a form that transforms naturally 
under the conformal group. 

The extensive part of the energy was computed in the previous subsection with the result 
(74). The volume factor in the static frame (y) is related to the infrared cutoff in the matrix 
model coordinates (x) through x ~ e y so we can express this result concisely in terms of the 
energy density 

e = -i < 76 > 

in the static frame. Note that this formula is exact at J + since the subleading terms in (73) do 
in fact vanish as Q + — > oo with Q_ fixed. The classical expression for the pressure density is 
identical to that of the energy density so, repeating the arguments in the previous subsection, 
we find p = e. The complete energy-momentum tensor in the static frame thus takes the form 

(Ty+,y+) = (Ty_, y _) = — ^g^> (77) 

in light-cone coordinates 

Let us now transform to the cosmological coordinates Q which contain the behavior at 
late times. Since (75) is a conformal transformation we can use the standard, anomalous, 
transformation rule 

Tq-,q- = (|^) 2 W + m {v ~' Q - }s ' (78) 

of the EM-tensor. The second term is present because the short distance singularity of the EM- 
tensor, although still cancelled by the explicit counter-term (65) in the collective field theory, 
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differs by finite amounts in the two coordinate systems. The expectation value of (78) relates 
the EM-tensor of the static and the cosmological coordinates. Evaluating the derivatives we 
find 



dy_ 1 

dQ- ~ 1 - \e-Q- 



(79) 



Ae"«-(1 - \\e 



<»-«->« = (1 -xA r ' m 

and then (78) becomes simply 

< W> = (81) 

This means physical observers will experience no outgoing flux of energy at X + . This result is 
a tremendous surprise. 

In standard computations, such as those considering moving mirrors, the static EM-tensor 
{Ty_, y _) would be taken to vanish, and so only the second term (the Schwarzian derivative) 
in (78) would contribute. The associated energy would be interpreted as the energy of the 
particles produced by the expansion. This is also the natural interpretation of the Schwarzian 
derivative here. However, in the present context the particle production is partially obscured 
by another effect which is operative as well: the (one-loop) energy density of the "in" vacuum 
is non-zero and finite in the co-moving frame (y). The corresponding EM-tensor is blueshifted 
in the frame of asymptotic observers (Q) on X + and this leads to an energy flux. This effect 
is the origin of the first term in (78). Our remarkable result is that these two effects cancel 
each other precisely such that the cosmological vacuum has precisely the same EM-tensor as 
the static vacuum. This is very different from the standard computations. 

The EM-tensor (81) takes the same value along the entire X+ as it does on Z~ . This suggests 
that, despite appearances, X + is a perfectly nice locus to define large classes of observables. One 
may consider the scattering of particles from X - to X + , and correlation functions on X + itself. 
The existence of such observables cannot be taken for granted in a time-dependent setting but, 
in matrix cosmology, it would seem that they both exist and are computable. It would clearly 
be interesting to make this more explicit. 

Our result that the one-loop energy of the matrix cosmology agrees with that of the static 
ground state is all the more surprising because the classical energy in this specific background 
does receive the very specific time-dependent contribution 

E cl = -F J ^ dxx^/x 2 - 1\i (82) 

where E^ is the static ground state energy. This is consistent with (33). That the one-loop 
energy is the same as in the ground state of course means there is no energy flux. 
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In summary, we have found a remarkable cancellation between two apparently different con- 
tributions to the EM-tensor. The consequence is that the outgoing EM-tensor T is identical 

in the in and out vacua. This is quite unusual and surprising, at least to us. It is natural 
to suspect that the cancellation ultimately stems from the integrable nature of the underlying 
matrix theory. To understand this better, it would be interesting to make the formal descrip- 
tion of the matrix cosmology as a deformed matrix model more concrete. It would also be 
interesting to extend the computation of the EM-tensor to the broader class of cosmological 
solutions generated in section 2.2 from the W^-symmetry of the model. 
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